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C^ ' We discover numerically that a moving wave packet in a quantum chaotic billiard will always 

^-H , evolve into a quantum state, whose density probability distribution is exponential. This exponential 

^3 . distribution is found to be universal for quantum chaotic systems with rigorous proof. In contrast, 

CN| ■ for the corresponding classical system, the distribution is Gaussian. We find that the quantum 

^.H \ exponential distribution can smoothly change to the classical Gaussian distribution with coarse 

^ , graining. 
^^ ■ 

(n: 

The study of quantum chaos intends to elucidate the intriguing correspondence between quantum mechanics and 

' ^ \ classical mechanics by focusing on the effects of chaos of a classical system on its quantum counterpart [l|-|5| . Intensive 

^ ■ efforts have found certain universal characteristics of quantum chaotic systems, for example, the Wimer distribution 

W)| of energy level spacing [2] and the "scarring" of eigenstates along the classical periodic orbits [6l, I?!]. By following 

■^^ ■ a pioneering work [8], we study the dynamical evolution of a moving wave packet inside a ripple quantum billiard 

^1 [y]. We find that the wave packet evolution always leads to an "equilibrium" state, where the density probability is 

^ . exponentially distributed. With rigorous proof, we show that this exponential distribution is universal for all quantum 

qH| chaotic systems. This exponential law is in stark contrast with the evolution of a "cloud" of classical particles, which 

^ . always leads to a Gaussian distribution. This universal behavior is observable with Bose-Einstein condensates (BECs) 

2 ! [lO^b-l] and other quantum systems. Furthermore, we demonstrate that the quantum-classical correspondence between 

Q . the exponential distribution and the Gaussian distribution can be achieved with coarse graining. 

_L I We consider the following Schrodinger equation with the units h = 2m = 1 

^^ ■ where V5 represents the hard- wall potential for a ripple billiard (see Fig. [B^O). The ripple billiards have the 

^ . advantage of having an exact Hamiltonian matrix in terms of elementary functions [9,]. The left and right walls 

^H ' are described by functions ^{h — acos (27r?//L)), respectively. The initial wave function is a moving Gaussian wave 

C""*^ . packet given by ^ {x^y^t = 0) = ^g (x^y) e^^o^/^, where ^g (x^y) is the ground-state wave function in a harmonic 

r^ ■ The time evolution of the wave packet density n (x, y^t) = |^ (x, ?/, t)\ is shown in Figs(Tl\l-A4 for the parameters 

^^ . ujx = ujy = 2^ vq = 10, L = 30, 6 = 15 and a = 6. The time is in the unit of T5 = 2 (a + b) /vq, which is the longest 
time for a classical particle making a round-trip along the x direction inside the billiard. As clearly seen in the figure, 
the wave packet moves and expands, and eventually gets reflected multiple times by the curved hard walls. As a 
result, the wave packet gets smeared out in the billiard. After about t = 6Tg, the density of the wave packet reaches 
kind of an "equilibrium" state, where the overall features of the wave packet no longer changes. For comparison, the 
^ ■ evolution of the same wave packet in a square billiard (a = 0) is shown in Figs(Tpl-B4. The difference is obvious: in 
the square billiard, the wave packet always has certain patterns, looking very regular, and it never settles into any 
fixed pattern. 

With the inspiration to describe quantitatively the contrasting images between Figs(T]A.3-A4 and Figs(lp3-B4, we 
introduce the density probability distribution function P (n, t), which is defined as the probability of the wave packet 
density being between n — Sn/2 and n + 5n/2. Its mathematical expression is 

o .^ .^ _ ^ (^ - "^V^, t)-S{n + 5n/2, t) 

on • ^total 

where 5'totai is the area of the whole billiard and S (n^t) is the area of the regions, where the density is larger than n. 
This density probability distribution is normalized by definition. As the wave packet in our calculation is normalized. 
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FIG. 1: Dynamical evolution of a moving wave packet in a ripple billiard and a square billiard. (A0-A4) The densities of the 
wave packet in the ripple billiard at different times t. (B1-B4) The densities of the wave packet in the square billiard at the 
same set of times except at t = 0. (C) The time evolution of the probability for the density being above the average. The 
red-dashed line is for e~^ with e being the Euler's number. (D) The density probability distribution functions at t = 12Ts. The 
red solid line is for exp(— no). The inset shows the semilog plot of the same sets of data. For the ripple billiard, the parameters 
are L = 30, b = 15, and a = 6. The initial speed of the wave packet is vo = 10. to is in unit of Ts = 2(a + b)/vo and no is in 
unit of the averaged density Us. h = 2m = 1 are adopted in the compuation. 6n = 0.02 is used to obtain the results in (D). 

the averaged density is n^ = l/S'totai- With the averaged density, we define a dimensionless density as no = n/us. 
The dimensionless probability distribution function is then Pq (no, to) = ngP in^t) with to = t/Tg. 



Figure [Tp shows the time evolution of Sait) = S {ris^t) /S'totab the probability of having density above the average. 
For the square billiard (a = 0), Sa{t) exhibits large-amplitude oscillations, reflecting the ever-changing regular patterns 
in the density. Very differently and also strikingly, Sa quickly reaches a plateau, fluctuating slightly around a fixed 
value, which is around e~^ with e ^ 2.718 being the Euler's number. This signals that the wave packet finally evolves 
into an "equilibrium" state, where the overall feature of the wave function no longer changes. The distribution 
function Pq (no, to) also has similar behavior: for the square billiard, Pq always changes with time; in contrast, for the 
ripple billiard, Pq settles into an "equilibrium" function. As shown in Fig. [Ip, we find by numerical fitting that the 
"equilibrium" function is exponential. 



Po^ino) 



(3) 



This simple exponential law is also found numerically for the the widely studied stadium billiard [6|, |8| . The evolution 
of a wave packet in Fig(T] can also be regarded as an evolution of a BEC without interaction. In fact, all the 
parameters chosen for Fig. [T]come from relevant BEC experiments [lO|. We find that for sufficiently large interaction, 
the density distribution of a BEC in the square billiard can also reach the above exponential distribution after 
long-time evolution. These numerical calculations indicate that the exponential law in Eq.(|3]) is likely universal and 
may apply in any quantum chaotic systems and in many different settings. To support this, we provide two more 
sets of numerical results. 






FIG. 2: Two more examples of exponential distributions. (A) Density distribution in the momentum space for the quantum 
state shown in Fig[T]A.4. (B) Probability distribution for the momentum distribution in (A). In obtaining (B), the momentum 
space is divided into concentric rings with equal spacing 5k = 7r/30 and each ring is then divided into equal pieces so that each 
piece has the area oi ^^ 5k x 5k. The density at each unit piece is renormalized with the averaged density of the ring where 
the unit piece belongs. The large fluctuations in (B) are caused by the small number of sample points. (C) The density of the 
randomly superposed plane waves. It is computed according to Eq.Q with 10,000 plane waves for an area of 400x400. (D) 
Density probability distributions for the random quantum state in (C). The distribution functions are obtained by dividing 
this whole area equally into different square cells, 0.1 x 0.1 (black line), 2 x 2 (green line), and 8 x 8 (red line). 



In FiglSl^, we have plotted the momentum distribution of the wave function shown in Fig(T]A.4. The momentum 
distribution has a ring-shaped structure with a spotty look. To understand this structure, let us consider the 
time evolution of a "cloud" of non-interacting identical classical particles, whose initial velocity distributions 
are identical to the wave packet in Fig(T]A.O, that is, G(l/,t = 0) ~ g-|^-^o| /2cr^ p^^ ^ chaotic billiard with 



hard- wall boundary (e.g., the ripple billiard), after long-time evolution, the velocity distribution will become 

G{v^t :^ Tg) ~ e v ' '^oj/ ^v^ which is exactly a ring-shaped distribution with a maximum at 'U = vq. However, 
G{vjt ^ Tg) is very smooth, implying that the spotty look in FiglSJA. has a quantum origin. As we have done for 
the density in the real space, we use similar statistics to quantify the spotty or random image in Fig|2]A.. We first 
divide the whole momentum space into a series of rings with equal small spacing and then cut each ring into small 
pieces of the same area. In this way, we obtain a set of momentum densities. This set of densities does not follow the 
exponential law. However, after renormalizing the densities of each ring with the averaged density for the same ring, 
we again find an exponential distribution as shown in FiglSB- 



In another aspect, following the spirit in Ref.[8|, we construct a semiclassical wave function with the classical 
momentum distribution G 



^ (x, y,t^Ts)oc f dv^dvyG (it, t > T,) e-^(^-^/2+^^ 



yy/2-\-ipriVa;,Vy,t)) 



(4) 



where (pr {vx^Vy^tys are assumed to be random phases caused by the classical chaos. Since G does not fully describe 
the quantum momentum distribution, the so-constructed wave function ^ should be different from the wave function 
shown in Fig. [T]A.4. This is indeed the case: the density computed with the above wave function (Figl2p) looks 
different in structure from Fig. [T]A.4 although both of them look random. However, despite this difference, the 



random wave function ^ in EqHJalso obeys the exponential law (black solid line in Fig|2p). Due to the large sampling 
size, the numerical result fits the exponential function exp(— no) almost perfectly. 

All the above numerical results point to the universality of the exponential distribution. This will become apparent 
as we give a rigorous proof in the following. We divide the region of the dynamical evolution into N equal parts. 
Then the quantum state |^ (t)) can be approximated as 

N 



Here aj{t) = J^ ^ {x^y^t) dxdy/^Siota\/N with Sj denoting the jth part. Obviously, Ofj's are complex numbers, 
satisfying the normalization condition 

El«.f = l. (6) 

We assume that for a fully chaotic classical system, its corresponding quantum dynamics will always drive the system 
to states, where the Ofj's are a set of random complex numbers that satisfy the above normalization condition. With 
these considerations, the probability of \aj\ being between 7^ and 7^ + ^7^ is 

P (7i) djj = ^ . ^ \ ^ ^d7i- (7) 

jd'^ai---d'^aNS\l-Y,i=i\(^i\ ) 

It is easy to find that J^ P (jj) djj = 1. After straightforward calculations |l2|, we have 

P(7,) = (iV-l)(l-7,f"'- (8) 



In the limit of TV ^ oo, we have aj{t) ~ ^(xj , yj , t)/ yJugN and the density probability distribution for the jth part 
becomes P {rij) = e~"^^, where rij = |^(xj, ?/j)p/ns. Because different part satisfies the same exponential probability 
distribution, the dimensionless density probability distribution Pq (no) defined for the whole region satisfies the 
exponential probability distribution e~^°. 

The above rigorous proof indicates that this exponential law is universal and should appear in any quantum chaotic 
systems and many different settings, for example, in the momentum space. The proof can be used straightforwardly 
or with small variation to understand three sets of numerical results presented before the proof. With the universality, 
it is not surprising to notice that this exponential law is known to exist in the speckle pattern of light [13]. However, it 
seems that this exponential law has been limited to this special system until now, and has never been fully discussed 
in a wide context or explored for broad applications. 

The established exponential law is in fact a quantum phenomenon. We can imagine to throw many, many classical 
particles into the billiard with equal probability to every part of the whole region. By the central limit theorem, 
the resulted density probability distribution is a Gaussian with its peak located at the averaged density, drastically 
different from the exponential distribution for a quantum gas. To confirm this, we have simulated the evolution of a 
cloud of non-interacting classical particles with the same spatial and velocity distributions as the wave packet studied 
in Figdl The results are shown in Fig|3l where we see the classical cloud expands and reflects very much like its 
quantum counterpart, the wave packet in Fig(TJ However, the density probability distribution function for the classical 
cloud after long time evolution is Gaussian (see Figl3p) as we just expected. As a result, we can say that the ex- 
ponential distribution is a quantum effect and it can be used in experiment to tell whether a gas is quantum or classical. 

Eigenstates of a quantum chaotic system have been widely studied. Many of these eigenstates tend to scar over 
classical periodic orbits [6]. This means that eigenstates have some regular pattern and it should not have the 
exponential distribution. We have checked this numerically and one example is shown in Fig HI In fact, the density 
probability distribution of an eigenstate was studied in literature in terms of amplitude probability distribution and 



no similar universal behavior was found [14 1. 
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FIG. 3: Evolution of a cloud of classical particles. (A0-A4) are plotted with 20000 particles while the statistics shown in (B) 
is done with 1 million particles. The cloud initially has the exact same spatial and velocity distribution as the initial wave 
packet in Fig[T] The unit cell used for obtaining the statistics in (B) is 0.6x0.6. 
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FIG. 4: The 880th eigenstate of the ripple billiard and its density probability distribution. The parameters for the ripple 
billiard are the same as in Fig[T] The inset in (B) is the semilog plot of the same data. 

The central question in the study of quan tum chaos is how the classicality observed in our daily life emerges 
from the underlying quantum world [l5|, ll6|. Therefore, it is very interesting to see how the classical Gaussian 
distribution emerges from the quantum exponential distribution. We find that the quantum-classical correspondence 
in this regard is built by coarse graining. We have computed the density of the random wave in Eq.(j4j) for an area of 
400 X 400. To obtain the density probability distribution Pq, we then divide the whole area equally into small square 
cells and sample the averaged densities in these small cells. We find that the distribution function changes its shape 
with the sizes of the cells: as the cell size increases, the distribution changes from exponential to Gaussian as shown 
in Fig|2)I). The same transition is observed for the density in Fig(T]A.4 with coarse graining although the fluctuations 
are large due to the small sampling size. Besides its fundamental significance, this coarse graining result has also 
implications in practice: with a camera of low resolution, one is likely to observe a Gaussian distribution even for a 
quantum gas. 



We have been calling the quantum state reached after long-time evolution "equilibrium" state. This is because we 
can not resist the temptation to compare it with the equilibrium state in thermodynamics. In thermodynamics, a 
system will always evolve into an equilibrium state after fully relaxed. After the equilibrium is reached, the state 



still changes microscopically; however, the corresponding macroscopic observables do not change, in particular, the 
statistical distribution function no longer changes. This is the same here for the long-time quantum state: after a 
certain long time (~ 6Tg in our case), the quantum state still changes locally while its overall structure remains the 
same and the density probability distribution function sticks to the exponential law. It is worthwhile to compare 
the hyper-sphere in the Hilbert space defined by the normalization condition (|6j) to the hyper-sphere in the phase 
space defined by a constant energy. As we know, for almost all the microscopic states in the phase space sphere, they 
share the same macroscopic characteristics, such as the same pressure and the same thermal expansion coefficient. 
Similarly, with the understanding that we have gained so far, almost all the states in the Hilbert space should have 
the exponential distribution. It is true that there are many states such as the eigenstate shown in FigHl which do not 
have the exponential distribution. However, we expect the measure of these states is zero. In light of this discussion, 
it is likely that this exponential distribution is related to the thermalization of isolated quantum systems [17]. 

Quantum chaotic system has been defined as the quantum system whose corresponding classical system is chaotic. 
With the exponential distribution law, we may be now ready to define quantum chaotic system without referring to 
classical system. One might choose to define a quantum chaotic system as a quantum system which has the ability 
to drive an initially regular wave packet to a coherent superposition of completely random wave functions, which has 
the observable exponential density probability distribution. We finally note that for the exponential distribution, the 
density fiuctuation is dm? = n^, which is much larger than the Gaussian distribution. In other words, the quantum 
fiuctuation in density is much larger than the classical fiuctuation. This large quantum fiuctuation in density may 
help to explain the formation of stars if a wave function for the whole universe is assumed [18.. .19.]. 
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